Interplay of roughness/modulation and curvature at proximity 
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We show that roughness or surface modulations change the distance dependence of (power-law) 
interactions between curved objects at proximity. The modified scaling law is then simply related 
to the order of the first non- vanishing coefficient of the Taylor expansion of the distribution of sep- 
arations between the surfaces. The latter can in principle be estimated by scanning measurements, 
or computed for well characterized modulations, and then used to predict short-distance scaling be- 
havior in disparate experiments. For example, we predict that the radiative heat transfer between 
a rough sphere and a plate approaches a constant with decreasing separation. Similar saturation is 
expected for the Casimir force between dielectric or metallic surfaces with appropriate modulations 
over distinct length scales. 

PACS numbers: 12.20.-m, 44.40.+a, 68.35.Ct 
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The Proximity Approximation (PA) has long served as 
a useful guide for estimating interactions between closely 
spaced objects with curved surfaces. Originally intro- 
duced by Derjaguin [l| to compute van der Waals forces 
between colloidal particles, the PA relates the interaction 
between curved objects at close separations to the corre- 
sponding interaction between two flat surfaces, over an 
area determined by the local radii of curvature. The PA 
has been successfully used in experimental studies in a 
variety of fields, such as quantum Casimir forces (H - !), 
classical Casimir forces in a fluid near a critical point 5[ , 
near field radiative heat transfer 0, 0] , and interactions 
among nuclei Q . Recently it has also been applied to the 
gravitational force in searches of non-Newtonian grav- 
ity 1- _ 

In this paper, we utilize PA to investigate scaling of 
the interaction between two surfaces whose local radii 
of curvature can be decomposed as the sum of compo- 
nents which vary along the surface on well separated 
length scales. We note that this situation encompasses 
the experimentally relevant cases of two large rough ob- 
jects, such that the correlation length of the roughness 
is much smaller than the characteristic radius of curva- 
ture of the surface, or alternatively the case of a sur- 
face with large average radius of curvature modulated 
by small structures fabricated by the experimenter. We 
find that at proximity a subtle interplay between rough- 
ness/modulation and the global curvature of the sur- 
faces leads to a drastic change in the distance depen- 
dence of the interaction, compared to that for perfectly 
smooth (structureless) surfaces. The modified scaling law 
is found to depend in a simple way on the order of the 
first non- vanishing term in the Taylor expansion of the 




FIG. 1. Profiles of two gently curved surfaces: The local sep- 
aration S is the sum of a slowly varying component S c due to 
overall shape, and a rapid part due to roughness/modulation; 
d is the distance of closest approach. 



height distribution function of the surfaces. The latter 
is a geometric feature of the surface that can be either 
computed for well characterized modulations, or experi- 
mentally obtained via scanning probes of rough surfaces, 
and then used to predict the short distance scaling of 
results in a multitude of experiments. We discuss the ex- 
ample of radiative heat transfer, and show how modula- 
tions (of various shapes) or roughness limit the maximum 
amount of non-contact heat transfer that can be achieved 
between two surfaces in close proximity. Applications to 
the quantum Casimir force shall also be discussed briefly. 

Consider two objects with surface-to-surface separa- 
tion S = S(x,y), where (x,y) are Cartesian coordinates 
in a plane that separates the two surfaces, as in Fig. [T] 
The unspecified interaction / between them may depend 
in a complicated way on the geometry of the objects. 
However, we assume that at close proximity it is dom- 
inated by the local properties of the surfaces near the 
point of closest approach. In the proximity approxima- 
tion Ipa the net interaction is given by 



Ipa 



dx \ dylpp (S(x,y)) , 



(1) 
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where I PP {S) denotes the interaction between two par- 
allel plates at separation S (per surface area). It is 
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reasonable to expect that Eq. ([I]) becomes exact for 
nearly flat surfaces, i.e. in the limit d/p — > 0, where 
p sets the scale for the local radii of curvatures. Indeed, 
in certain cases Eq. (TTJ represents the leading term of 
an expansion of the interaction in powers of the gra- 
dient of S [lMl (see Eq. CLll below). This implies 
that Eq. ([1]) is indeed asymptotically exact in the limit 
d/p — > 0, the leading curvature correction being of higher 
order in d/p. Equation ([I]) can be alternatively re- 
cast in terms of the height distribution function fLU Il5j 
f(s) = J dxdy 6(S(x,y) — d — s), as 

Ipa = J dsf(s-d)I pp (s). (2) 

Note that f(s) is measured with respect to the distance 
d of closest approach, such that /(s) = for s < 0. It 
is thus a function of shape (and orientation) of the ob- 
jects, and independent of their separation. Equation ([5]) 
highlights the basic assumption of PA in neglecting the 
detailed topology of the surface (e.g. gradients), being 
sensitive only to the projected surface area at separation 
s. While in general only an approximation, Eq. (|2|) can 
be used to understand the qualitative effects of modu- 
lation or roughness in concrete experimental setups; the 
function /(s) can be easily extracted from atomic force 
microscope (AFM) scans of the surfaces [l4j|, or com- 
puted analytically for model surfaces (e.g., below). 

In what follows, we assume that I PP (d) diverges with 
a simple power-law for cZ — 5- 0, 

Ivp(d) = ^. (3) 

For example, v = 3 for the quantum Casimir energy be- 
tween perfectly conducting surfaces Q, while v = 2 for 
the critical Casimir energy involving binary fluids [j| or 
radiative heat transfer [16j |. The dimensional coefficient 
a depends on the specific system under study. 

Starting from Eq. ([3]), it is easy to verify that the 
asymptotic behavior of Ipa at close proximity is deter- 
mined by the first non- vanishing coefficient f( n \ff) of the 
Taylor expansion of / near s = 0, as summarized in Ta- 
ble HI Considering for example a sphere of radius R in 
front a plate, we have 

f s (s) = 2v(R-8), (4) 

for < s < R and zero otherwise; /s(s = 0) is finite 
(case C = 1 in Table [TJ) , leading to well known results: 
Ipa oc d~ 2 for v = 3, and Ipa oc d^ 1 for v = 2. Another 
example is a square base pyramid of height h and base 
length I, with axis perpendicular to a plate, for which 

f P (s) = 2 S i,, (5) 

for < s < h and zero else. The sharp tip of the pyra- 
mid has vanishing surface density (case 2), such that PA 
yields a distinctly different behavior than for the sphere- 
plate configuration: A I/d-divcrgcnce for v = 3, or a log- 
arithmic one for v — 2. Although the PA is not asymp- 
totically exact for singular surfaces with sharp tips, the 
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TABLE I. Asymptotic behavior of the interaction between 
two surfaces, as determined by the order of the first non- 
vanishing Taylor coefficient /'*^ (0) of the height distribution 
function of the surfaces. The interaction scales as I PP (d) = 
a/d v for two parallel surfaces. The length do depends on 
higher derivatives of /. 

above conclusion agrees qualitatively with exact compu- 
tations of the Casimir force and heat transfer in the cone- 
plate geometry [It], HH • 

We now turn to the experimentally relevant case of 
rough or modulated surfaces. The surface-to-surface sep- 
aration S can be decomposed as S = S c + S r , where S c 
is the separation between two smoothened surfaces with 
general radius of curvature 1Z, and S r is the additional 
shift resulting from roughness or modulation (see Fig. [I]) . 
If the length scale £ characterizing the roughness or the 
modulation is much smaller than the characteristic size L 
of the interaction area (typically of order L ~ s/dlZ), i.e., 
£ <C L <C 7Z, we find that the height distribution func- 
tion f(s) for S(x, y) is well approximated by the convolu- 
tion of the corresponding functions f c of the smoothened 
surface, and f r (normalized by a unit surface area) for 
roughness/modulation, i.e. 

f(s)= f d S 'f c (s')f r (s-s'). (6) 
Jo 

This result also obtains from the commonly-used two- 
step application of the PA in which the PA is first 
used to estimate the correction arising from roughness 
or modulation (per unit area) to the interaction between 
parallel plates, and then used again to estimate the ef- 
fect of curvature. It has been noted (see e.g. Ref. (f5j ) 
that if the experimental conditions are such that the ap- 
proximate formula of Eq. © is not valid for the surfaces 
of interest, measurements might not be fully repeatable, 
as they would be sensitive to e.g. uncontrollable relative 
displacements of the two surfaces in the (x, y) plane. 
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The asymptotic scalings in Tabic U also apply to the 
combined distribution function in Eq. ([5]) , where we note 
an intriguing relation. Restricting to regular, analytic 
functions near s = 0, the case number of the combined 
function / in Eq. (|6]) is precisely the sum of the case 
numbers of the individual functions f c and f r , 



C(f)=C(f c )+C(f r ). 



(7) 



Equation ([7]) can be derived with the convolution theo- 
rem for Laplace-transforms, and combined with Table U 
provides a means for obtaining the modified scalings. 

Figure[2]shows f(s) for the case of a sphere with differ- 
ent types of modulation in front a plate where the 
red curve is for the smooth sphere of Eq. ((4]). If the sur- 
face of the sphere is covered with (square base) domes of 
height h (f r = fu = 2(h — s)/h 2 ), f is found analytically 
with the help of Eqs. Q and (JBJ, as 



f SD (s < h) = 



2ns (6hR 



3hs-3Rs + s 2 ) 

-w (8) 



As we see, /sd(0) vanishes, however its first derivative 
does not, see the blue curve in Fig. [2] Therefore, 
the dome-like modulation (C = 1), yields C = 2 when 
convoluted with the sphere-plate geometry. 

Suppose now that the sphere is covered with pyramids 
of height h (again square based for ease of tiling, as in 
Eq. ©). After convolution with the sphere we now find 



fsp(s <h) = 



2irs 2 (3i2 - s) 
Zh 2 



(9) 



which is order s 2 (green curve in Fig. [2]). Sharp tips (C = 
2), give C = 3 when convoluted with a curved surface. 

We now consider surface roughness, which in general 
depends on materials and manufacturing (20l - [22l | . While 
small deformation of amplitude w -C d can be treated 
perturbatively [23],|24[, perturbative methods fail on close 
approach when d becomes comparable to w, and the PA 
provides a better route [25| . For practical ease, the height 
distribution function f r of a rough surface will be approx- 
imated by a Gaussian [TJ[, 



fn(s > 0) 



A/"crV27r 



exp 



(g - go) 
2a 2 



(10) 



where a is the standard deviation and Af is chosen such 
that J °° dsfn(s) = 1. The distribution f r has no sharp 
boundaries (as the peaks and grooves of the surface can in 
principle have arbitrary extension) , and we introduce an 
additional parameter sq to set s = 0, i.e., the "touching 
distance." This parameter sets the height of the highest 
peak, which is touched first on approach and acts as a 
shift on the separation axis. The truncated Gaussian in 
Eq. (p~0|) corresponds to case 1, although with a possibly 
small value of /i?(0). After convolution with the curva- 
ture of the sphere, we find that the height distribution 
function of the rough sphere is linear near s — (C = 2). 




FIG. 2. Height distribution functions for model surfaces (in- 
dicated by icons/text) in the sphere-plate geometry; for mod- 
ulated surfaces h — R/10, while a — i?/40 and so = 2a for the 
rough surface. The larger range of s in the inset emphasizes 
that the impact of deformations is mostly for small s. 




d [nm] 

FIG. 3. Heat transfer between a flat plate and deformed 
model spheres (as shown by icons/text) of radius R=50/im. 
The amplitude h for modulations (dome or pyramid), and the 
width a for roughness are as indicated, and 4 times smaller in 
the inset. The two curves for roughness correspond to so = 2a 
(upper) and so = 3a (lower). 



The magenta curve in Fig. [5] starts linearly, but with a 
small slope (for the chosen value of so = 2a), such that 
the rough surface is intermediate between the dome-like 
and pyramid-like modulations considered earlier. 

As a specific example, we consider radiative heat trans- 
fer H [26[ between a rough or modulated sphere and a 
(perfectly flat) plate, demonstrating that roughness or 
modulation limit the maximum achievable non-contact 
energy transfer, an issue of much current technological 
and experimental interest 0, 0, H3 ■ We note that for 
heat transfer, the PA has been shown to approach the 
exact result at small separation for a sp here in front a 
plate \2§l . [29| . and for two spheres |30j . As shown in 
Ref. [l3j, an asymptotic expansion in powers of the gra- 
dient of S exists for heat transfer H, which for small 
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FIG. 4. Example of a surface pattern (not to scale) with 
C — 4. For such a geometry, the quantum Casimir force for 
real materials (and in general any interaction with v = 3) is 
predicted to level off to a constant value at small separations. 

separations leads to 

H = H PA + j ds p g(s - d)H pp (s) + ■■■ . (11) 

The second term in Eq. (jlip represents the leading cur- 
vature correction to PA, with g(s) = J dxdy S(S(x,y) — 
d — s)VS ■ VS 1 encoding the mean squared gradient at s. 
The dimensionless coefficient (3 depends on material and 
temperature, and is typically of order unity The 
function g(s) for a rough or modulated surface can be 
found in a manner similar to Eq. (JSJ), and then Eq. (|11[) 
allows us to judge whether the scaling laws obtained from 
the PA are asymptotically exact for specific shapes. 

For small separations d, the radiative transfer between 
parallel plates (per unit area) diverges as a power law 
with exponent v = 2 [3l[ • Figure [3] illustrates the trans- 
fer as given by PA for the cases shown in Fig. [2] In 
all cases we subtracted from H the far field contribution 
(computed at the reference separation d = 300 ran) . The 
smooth sphere shows the expected divergence of 1 jd (red 
curve); dome-like modulations (with /(s) in Eq. ([5])) re- 
duce this divergence into a logarithmic form (blue curve), 
such that 

R d 

lim H SD (d) = -Ana- log — . (12) 
d->o h do 

For this type of modulation, the gradient correction in 
Eq. (ITT1) is of order d° , and therefore Eq. (IT2^) becomes 
exact for vanishing separation. 

A surface with pyramid-like modulations (/(s) as per 
Eq. (|9])) shows an even more drastic slow down: accord- 
ing to Table HI the transfer levels off to a constant value 
(green curve). However, evaluation of the gradient term 



in Eq. (|lip shows that this correction is now of the same 
order as the PA term. From this we infer that PA is 
not exact in this regime, while the predicted saturation 
to a constant is presumably correct. For rough surfaces 
the height distribution function is somehow intermediate 
between those for dome-like and pyramid-like modula- 
tions. Confirming this expectation, we find that H di- 
verges logarithmically for rough surfaces, but only with 
a small prcfactor. Depending on a and sq in Eq. (|10l) . 
this divergence may not be visible in practical cases, 
giving the appearance of saturation as in the case of 
pyramid-like modulations. Reducing the amplitude of 
roughness/modulation does not change the asymptotic 
behavior, but differences from a smooth sphere become 
visible only at closer separations (inset of Fig. [3]). 

On a practical note, the right axis in Fig. [3] shows the 
ratio to the classical far field (d — > oo) transfer for SiC>2, 
for Ti = and T 2 = 300 K (which is - 4.2/iW [II]). 
With a = 0.2558nW [Hj], and taking R = 50/im and 
ti = 10 nm, we predict that (rather disappointingly) heat 
transfer can only be approximately doubled by reducing 
the separation to 1 nm (for the considered range of sq). 

In principle saturation can also be achieved for Casimir 
forces with v > 3 by appropriate surface fabrication. 
Figure @] depicts a combination of three height profiles, 
a sphere with a smooth modulation and an additional 
sharp modulation. If the corresponding three lengths are 
well separated, e.g. for R = 100/im, and modulations 
lengths of 1/im and 10 nm, respectively, the convolution 
in Eq. (JSJ) (and also Eq. j?}) can be applied two suc- 
cessive times, to yield C = 1 + 1 + 2. The Casimir force 
between two such dielectric or metallic objects, for which 
v = 3 3, should saturate to a constant value. 

We have shown that the interplay between roughness 
or modulation and curvature leads to drastic modifica- 
tions of the short distance behavior of interactions be- 
tween two surfaces. By simple computations based on 
the PA, we find that the scaling law for the interaction 
is determined by the first non-vanishing Taylor coeffi- 
cient of the height distribution function of the surfaces. 
If modulations occur over widely distinct length scales, 
their effect on the predicted scaling law is additive. 
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